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Abstract: This paper studies exact linearization methods for stochas¬ 
tic SISO affine controlled dynamical systems. The systems are defined 
as vectorfield triplets in Euclidean space. The goal is to find, for a 
given nonlinear stochastic system, a combination of invertible trans¬ 
formations which transform the system into a controllable linear form. 
Of course, for most nonlinear systems such transformation does not 
exist. 

We are focused on linearization by state coordinate transformation 
combined with feedback. The difference between Ito and Stratonovich 
systems is emphasized. Moreover, we define three types of linearity of 
stochastic systems — (^-linearity, cr-linearity, and gcr-linearity. 

Six variants of the stochastic exact linearization problem are studied. 
The most useful problem — the Ito - ga linearization is solved using 
the correcting term, which proved to be a very useful tool for Ito sys¬ 
tems. The results are illustrated on a numerical example solved with 
help of symbolic algebra. 

Keywords: exact linearization, feedback linearization, nonlinear dy¬ 
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1 Introduction 


The theory of exact linearization of deterministic dynamical systems has 
been thoroughly studied since seventies. This paper attempts to apply some 
of the results to the stochastic area. We emphasize the exact linearization by 
state coordinate transformation combined with feedback (further abbrevi¬ 
ated as SFB linearization). Our main goal is to identify the main difficulties 
of this approach and to consider applicability of the methods known from 
the deterministic systems. 

The task of SFB linearization is following: given a dynamical systems 0 
we are looking for a combination of coordinate transformation Tt and feed¬ 
back which will make the resulting system J~a,i3 ° "^(0) linear and 

controllable. One can also define the feedback-less linearization by coordi¬ 
nate transformation only (here abbreviated as SCT) or several variants of 
the input-output linearization. These variants are not considered here. 

Remark 1. The notation for composition of mappings sometimes differs; 
right-to-left convention is used here: / o g(x) := f{g{x)). V 

The subject of exact linearization of stochastic controlled dynamical sys¬ 
tems lies on the intersection of three branches of science: differential geom¬ 
etry, control theory, and the theory of stochastic processes. Each of them is 
very broad and it is virtually impossible to cover all details of their combina¬ 
tion. Hence it is necessary to choose a minimalistic simplified model for our 
problem and to refrain from most technical details. We decided to represent 
dynamical systems under investigation by triplets of smooth vectorfields and 
to concentrate on transformation rules for these triplets. The detailed inter¬ 
pretation of the vectorfield systems (i.e. solvability of underlying differential 
equations, properties of flows and trajectories) will be considered only on 
an informal, motivational level. 

For simplicity, we shall confine all the definitions of geometrical object 
to the Euclidean space; we will work in a fixed coordinate system using 
explicit local coordinates, which may be considered to be local coordinates 
of some manifold. This is mainly because we are unable to capture all 
consequences of the modern, coordinate-free differential geometry to the 
stochastic calculus (see eg. Kendall |^, Malliavin |]l^, Emery [^). We 
believe that this approach is quite satisfactory for the majority of practical 
applications. 
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1.1 Dynamical systems 

Definition 2. In this paper, a stochastic dynamical system 0 := (/, g, a, U, xq) 
is defined to be a triplet of smooth and bounded vectorfields /, g, and a de¬ 
fined on an open neighborhood 1/ of a point xq G M”’. We usually call U G MT" 
the state space, / the drift vectorfield, g the control vectorfield, and 
a the dispersion vectorfield. ■ 


From now on, let’s assume that all functions, vectorfields, forms, and 
distributions are smooth and bounded on U. 

In this paper, we will study almost only SISO systems, but in the case of 
stochastic MIMO systems with m control inputs and /c-dimensional noise the 
symbols g and a stand for nxm (nxk respectively) matrix of smooth vector- 
fields having its rank equal to m {k respectively). The class of all determinis¬ 
tic n-dimensional dynamical systems with m inputs will be called Hjyin, m) 
and the class of stochastic systems with fc-dimensional noise will be denoted 
with X(n, m, k). 

Similarly, autonomous deterministic dynamical system corresponds to 
a single vectorfield and a controlled deterministic dynamical system corre¬ 
sponds to a vectorfield pair. It is obvious that this approach is limited to 
time invariant, affine systems. 

Remark 3. The acronyms SISO and MIMO are used in the usual mean¬ 
ing even for systems without outputs, where the wording “scalar-input”and 
“vector-input” will be appropriate. Stochastic systems with m = 1 and k = 
1 will be considered SISO. V 


The definition may be interpreted as follows: there is a stochastic pro¬ 
cess xt defined on M"" which is a strong solution of the stochastic differen¬ 
tial equation dxt = f{xt)dt + g{xt)u(t) dt + a{xt)dwt, with initial condi¬ 
tion Xq, where u{t) is a smooth function with bounded derivatives and wt is 
a one-dimensional Brownian motion. The differential dwt is just a notational 
shortcut for the stochastic integral. 

Details of the theory of stochastic processes are beyond the scope of this 
article. The reader is referred to Wong and Hajek |Q, Oksendal |Q, Sage 
and Melsa |17|, Malliavin [12|, Kendall |^], Karatzas and Shreve Q; the text 
of Kohlmann [^] is freely available on the Internet. 

Theory of stochastic processes offers several alternative definitions of 
the stochastic integral, among them the Ito integral and the Stratonovich 
integral; each of them is used to model different physical problems. Con¬ 
sequently there are two classes of differential equations and two alterna¬ 
tive definitions of a stochastic dynamical system — Ito dynamical systems 
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defined by Ito integrals and Stratonovich systems defined by Stratonovich 
integrals. 

Serious differences between these integrals exists but from out point of 
view there is a single important one: the rules for coordinate transformations 
of dynamical systems defined by ltd stochastic integral are quite different 
from the transformation rules which are valid for Stratonovich systems. 

The definition of the Ito dynamical system used by us is formally equiv¬ 
alent to the definition of the Stratonovich system; the only difference will 
be in the corresponding coordinate transformation. 

If necessary, Ito and Stratonovich dynamical systems will be distin¬ 
guished by a subscript: 0 / G Xj{n,m, k) and ©5 G Xs{n,m, k). 

Remark 4. In this paper we will use the adjectives Ro and Stratonovich 
rather freely. For example we will speak of ’Stratonovich linearization’ 
instead of ‘exact linearization of stochastic dynamical system defined by 
Stratonovich integral’. V 

1.2 Transformations 

Furthermore, we will study two transformations of dynamical systems: the 
coordinate transformation Tt and the feedback The definition of these 

transformation should be in accord with their common interpretation. This 
can be illustrated on the definition of the coordinate transformation of 
a deterministic dynamical system Tt : X£)(n,m) —> X£)(n,m) which is 
induced by a diffeomorphism T: U ^ M"" between two coordinate systems 
on an open set U C M”. The mapping Tt is defined by: 

TT{f,g,U,xo) := {T,f,Tg,T{U),T{xo)). (1) 

Recall that the symbol T* stands for the contravariant transformation {T^,f)i 
Moreover, we will require that the coordinate transformation T 
preserves the equilibrium state of the system i.e. T{xq) = 0. 

The definition captures the contravariant transformation rules for differ¬ 
ential equations known from the basic calculus. 

Note that the words “coordinate transformation” are used in two differ¬ 
ent meanings: first as the diffeomorphism T: U ^ between coordinates; 
second as the mapping between systems Tt : Xo(n,m) ^ X£)(n, m). 

Coordinate transformation of stochastic systems distinguish between Ito 
and Stratonovich systems. One of the major complications of the lineariza¬ 
tion problems for Ito systems is the second-order term in the transformation 
rules for Ito systems: 
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Figure 1: Regular State Feedback 

Definition 5. Let U G M”" be an open set and let T: ^ M” be a diffeomor- 

ptiism from U to M” with bounded first derivative on U such that T{xq) = 0. 
The mapping Tt'- ^i{n,m,k) — > ^i{n,m,k) will be called a coordinate 
transformation of an Ito dynamical system induced by diffeomor- 
phism T if the systems 0i := {f,g,a,U,xo) and 02 := ^/, g,a,T{U),xoy, 
02 = Tr (0i) are related by: 

/ = TJ + PaT (2) 

m = T*gi for 1 < i < m (3) 

di = for 1 <i <k. (4) 


The symbol P^T stands for the Ito term which is a second order linear 
operator defined by the following relation for the m-th component of PaT, 

1 < m < n 

PaTm := 2 E (5) 

i J=1 * l=l 

The transformation rules for Stratonovich systems Pp'. X 5 (n, m, fc) —> 
Xs{n,m,k), {f,g,a,U,XQ) {Pf,Pg,Pa,T{U),T{xo)) are equivalent to 
rules valid for the deterministic systems; only the rule(^) for the drift vec- 
torfield must be added. 

The difference between the coordinate transformation of Ito and Stratonovich 
systems should be emphasized: in the Stratonovich case all the vector- 
fields transform contravariantly; on the other hand, in the Ito case, the 
Ito term PfjT is added to the drift vectorfield of the resulting system. 

Another important transformation of dynamical systems is the regular feed¬ 
back transformation. A feedback transformation is determined by two smooth 
nonlinear functions a: M” —and (3 : M"" ^ x M”* defined on U with (5 
nonsingular for every x £ U (see Figure |^. Usually, a is written as a col¬ 
umn m X 1 matrix; /3 as a square m x m matrix. 
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Definition 6. Let 0 = {f,g,a,U,xo) G X(n, m, A:) be a stochastic dynami¬ 
cal system. A regular state feedback is the transformation J-a,i3 ■ X(n, m, k) 
X{n,m,k) , {f,g,a,U,xo) ^ {f + ga,gP,a,U,xo). ■ 

A new input variable v is introduced by the relation u = a + Pv. Given 
the feedback J-a,i3 with nonsingular /3, we can always construct an inverse 
relation Ta,b •= ^ such that J^a,i3 ° ^a,b = ^a,b o ^a,i3 is the identity. The 

coefficients are related as follows: /3 = b~^, a = —b~^a, and a = —P~^a. 

This definition of feedback transformation can be used also for deter¬ 
ministic systems provided that the drift vectorfield a is assumed to be zero. 

The symbol JT,a,i3 is used to indicate the combination of coordinate 
transformation with feedback JT,a,i3 '■= ° can be interchanged. 

Remark 7. Observe that the order of feedback and coordinate transforma¬ 
tion in the composed transformation JTT,a,/3 '■= ^a,f3 ° '^T 

JT,a,f3 = Tto {f , g, U, xq) = Tt {f + ga, gP, U, xq) 

= {T*f + T^ga , T^gP, U, zo) = {T^f + {T^g)a , {T^g)P ,U, zq) 

= kFc,' o Tt- (6) 

The functions a(z)P p{z)' are equal to a{x) and P{x) written in the z co¬ 
ordinates a(z)' = a{x) o T~^{z), P{z)' = P{x) o T~^{z). V 

1.3 Linearity 

The definition of linearity is straightforward in the deterministic case. In 
contrast, the stochastic case is more complex, because there are two “input” 
vectorfields and thereby several degrees of linearity can be specified. 

Definition 8. The deterministic dynamical system 0£) = (f,g,U,0) G 
Xo{n,m) is linear if the vectorfield / is a linear mapping without no con¬ 
stant term and the vectorfields gi are constant; that is they can be written 
as f{x) = Ax, g{x) = B with A a square n x n matrix and B an n x m 
matrix. The matrices must be constant on whole U. ■ 

Definition 9. The stochastic dynamical system 0 = (/, g, a, U, 0) is: 

• ^(-linear if the mapping f{x) = Ax is linear without constant term 
and g{x) = B is constant on U. 

• cr-linear if the mapping f{x) = Ax is linear without constant term 
and (j{x) = S is constant on U. 
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• ^(T-linear if it is both g'-linear and cr-linear. 

The matrices A and B have the same dimensions as in Definition S is 
an n X k matrix. ■ 

Remark 10. The vectorfield g is constant on U if the value of g{x) is the 
same for every x & U. The vectorheld / on [/ is linear withont constant 
term if /(xq) = 0 and the superposition principle /(X 1 +X 2 ) = /(xi) + /(x 2 ) 
holds for every xi, X 2 ^ U. V 

We study systems at equilibrium i.e. we require that /(xq) = 0 and 
that all transformations preserve the equilibrium: T(xo) = 0, a(xo) = 0, 
and /3(xo) is nonsingular. The Ito systems require an additional condi¬ 
tion /(xo) + corra-(xo) = 0. The non-equilibrium case can be easily handled 
by extending the linear model with a constant term. 

Moreover we require that the resulting linear systems are controllable. 

A controlled deterministic dynamical linear system ©d = (Ax, i3, M”", 0) G 
X£)(n, m) is controllable if its first n repeated brackets form an n-dimensional 
space 


dim 


I 0<A:<n—l| = n. 


(7) 


Other definitions of controllability of linear systems exist. For example The¬ 


orem 3.1 of Zhou et al. [^| gives six definitions with proofs of equivalence. 

The controllability property deserves some attention in the stochastic 
case. The linear stochastic dynamical system is characterized by two input 
vectorfields g{x) = B and ct(x) = S. 


1. The definition of the controllability for the control vectorfield g is 
identical to the deterministic case; i.e. must be satisfied. This 
property will be called ^-controllability. 

2. We will also define u-controllability as the requirement that the 
repeated brackets S, AS, A^S ,..., A^~^S form an n-dimensional space. 

3. Finally, the linear system is gur-controllable if 

dim A*^i?, 0<A:<n — l||^ |a^5, 0 < k < n — = n. (8) 


In this paper, we do not deal with the reachability, controllability, accessi¬ 
bility, observability, and similar properties of nonlinear systems. 
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Figure 2: Dynamical System with a Prefilter. 


Definition 11. Let 0 = {f,g,a,U,xo) € X(n,m,A:) be a dynamical system 
such that /(xo) = 0. We call the combination of a coordinate transforma¬ 
tion Tt and a regular feedback Ta,i3 such that T{xo) = 0, q;(xo) = 0, and 
/3(xo) is nonsingular the linearizing transformation of 0 if the transfor¬ 
mation ° Tt converts 0 into a controllable linear system. 

For stochastic system we distinguish: 

(i) ^-linearizing transformation which transforms 0 into a g'-linear 
and ^-controllable system 

(ii) u-linearizing transformation which transforms 0 into a fi-linear 
and cr-controllable system 

(hi) ^(T-linearizing transformation which transforms 0 into a ^rcr-linear 
and ^-controllable system 

Note that for gicr-linearization we require g(-controllability. This is 
slightly stricter requirement than ga controllability but it should be 
naturally fulfilled by the majority of practical control systems. This 
requirement cancels many “uncomfortable” linear forms. Consider for 
example the system with prefilter of Figure |2| which is guT-controllable 
but g(-uncontrollable. 


The system 0 is linearizable if there exists linearizing transformation of 0. 
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1.4 Computational Issues 

In most practical circumstances, computational issues are the limiting factor 
of any application of differential geometric methods in control. 

The equations of exact linearization algorithms must be dealt in a sym¬ 
bolic form. Even the simplest exact linearization problems are extremely 
complex from the computational point of view. Therefore, the computer 
algebra tools are often employed. The results presented in this paper were 
tested by the author on few simulations of control systems in the symbolic 
system Mathematica . 

Of course, the computer algebra has apparently serious limitations and 
drawbacks. Viability of the symbolic computational approach to the prob¬ 
lems of the nonlinear control is studied by Jager p. Some very useful 
theoretical notes about the symbolic computation can be found in Winkler 
[p^ . Unfortunately, the limited scope of this article does not allow deeper 
discussion of these subjects. 

1.5 Applicatious 

We propose, very briefly, two applications of the theory presented here: 

(i) Control — a dynamical systems 0 obtained by exact linearization 
will be controlled using the linear feedback law: 

V = Kz -|- Kv, (9) 

where A is a row matrix of feedback gains, k is an input gain, z is 
the state vector, v is the original control input, and zz is a new control 
input. 

Two approaches can be studied — classical linear control methods and 
the more sophisticated stochastic optimal control approach studied for 
example by 0ksendal 0- 

The gKT-linear systems are natural candidates for such approach be¬ 
cause the other linear forms leave certain part of the resulting system 
nonlinear. 

(ii) Filtering — the hltering problem is probably the most useful appli¬ 
cation of the theory of stochastic processes. We want to give the best 
estimate to the state of a dynamical system defined by the stochastic 
differential equation: 

dxt = f{xt)dt + Of{xt)dwtf] (10) 
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based on observations of the from: 

yt = h{xt) + ah{xt) dwth. (11) 

xt is an n-dimensional stochastic process, /, u/, Uh are smooth vec- 
torfields and h is a smooth function. 

It would be interesting to use exact linearization of the nonlinear sys¬ 
tem to design an exact linear filter. Unfortunately, this idea has no 
direct association with the linearization results presented below, be¬ 
cause it requires output exact linearization or linearization of an au¬ 
tonomous system. Therefore, it would be helpful to extend our results 
to these cases in the future. 


1.6 Previous Work 


The problem of SFB ^-linearization of SISO dynamical system defined in 


the Ito formalism has been studied by Lahdhiri and Alouani []T|. The au¬ 
thors derive equations corresponding to (|108|) , ( 109 ) (eq 14, 15, 16 in |Tl|). 
These equations are combined and then reduced to a set of PDEs of a sin¬ 
gle unknown function Ti. Because there is no commuting relation simi¬ 
lar to (|^ the equations contain partial derivatives of Ti up to the 2n- 
th order (eq 23 in ||ll|| ). Next, the authors propose a lemma (Lemma 1) 
that identifies the linearity conditions with non-singularity and involutiv- 
ness of {adjgf, 0 < i < n — 2}. Unfortunately, we disagree with this result. 

It can be easily verified that for cr = 0 this statement does not correspond 
to the deterministic conditions (Proposition 0) , because the deterministic 
case requires non-singularity up to the (n — I)-th bracket, not only up to 
the (n — 2)-th one. Second, although the method of finding Ti was given 
(solving PDE), we do not think that the existence of Ti was proved. 

After this paper was finished, we discovered recent works of Pan [^] 
and 1^]. In the article |^] Pan defines and solves the problem of feedback 
complete linearization of stochastic nonlinear systems. In our terminology, 
this problem is equivalent to SFB MIMO input-output Ito ga linearization 
which was not studied by us. 


In I IT ] Pan declares and proves so called invariance under transformation 


rule which is exactly equivalent to our Theorem 18 which is probably the 
most important result of our paper. 

Althougth the problems solved by Pan were slightly different, he uses the 
same equivalence — Theorem This proves that our conclusions about 
applicability of the correcting term are perfectly valid. 
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In Pan consider three other canonical forms of stochastic nonlin¬ 
ear systems, namely the noise-prone strict feedback form, zero dynamics 
canonical form and observer canonical form also not studied by us. 


2 Deterministic Case 


In this section we recapitulate the results of the SFB and SCT exact lin¬ 
earization theory for SISO systems. For detailed treatment and proofs we 
refer to existing literature, above all the classical monographs of Isidori |^] 
and Nijmeijer and van der Schaft [O. For a very readable introduction to 


the field we refer to the seventh chapter of Vidyasagar |18|. The books also 
contain extensive bibliography. The monograph of Isidori Q builds mainly 
on the concept of relative degree. In contrast we will emphasize the approach 
of Vidyasagar [^] because the method is more suitable for the stochastic 


case. 


2.1 Useful Relations 

The solution of the SFB linearization problem as presented here uses the 
Leibniz rule 

with f,g smooth vectorfields on U] a: U —> M is a smooth function. The 
recursive form of the Leibniz rule allows to simplify the chains of differential 
equations for the transformation T. This can be expressed in the form of 
the following statement; 

For all X € U, k >0 these two sets of conditions are equivalent: 


(i) Cga = CgCfa = ■ ■ ■ = CgC'ja = 0 (13) 

(ii) Cga = ga = ■■■ = = 0. (14) 


Recall that the symbol Cjh stands for the Lie derivative defined hy Cfh = 
(/,grad/i) = Y17=ifiix) Higher order Lie derivatives can be defined 

recursively C^jh = h, = CjC^h for /c > 0. The Lie Bracket is defined 

as [f,g] := — ^g] there is also a recursive definition: 


adjg := g] 


ad^+^S- := 


/, ady g 


for A: > 0. 


(15) 
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Another very important result of the differential geometry is invariance 
of the Lie bracket under the tangent transformation T* (see Nijmeijer and 
van der Schaft Proposition 2.30 p. 50); 

Let T: U ^ M” be a diffeomorphic coordinate transformation, and / 
and g be smooth vectorfields. Then 

T4f,g] = [nf,ng]. (16) 


2.2 SFB Linearization 


Every controllable linear system may be, by a linear coordinate transfor¬ 
mation, transformed to the controllable canonical form (Kalman et al. |Q]). 
Furthermore, this controllable canonical form can be always transformed 
into the integrator chain by a linear regular feedback. Therefore, the inte¬ 
grator chain is a canonical form for all feedback linearizable systems. (See 
Vidyasagar []^| section 7.4). Consequently, the equations of the integrator 
chain can be compared with the equations of the nonlinear systems and the 
following proposition can be proved: 


Proposition 12. There is a SFB linearizing transformation ffT,a,i3 of a 
SISO deterministic dynamical system Qd = {f,g,U,xo) G X£)(n, 1) into a 
controllable linear system if and only if there is a solution Ti, T 2 ,..., Tn : M” —> M 
to the set of partial differential equations defined on U 


CfTi = Tj+i for 1 < i < n - 1 

(17) 

CgTi = 0 for 1 < i < n — 1 

(18) 

C-gTn 0. 

(19) 

Then the feedback is defined as follows: 


^f'^n n 1 

CgTn ^~CgTn 

(20) 

Proof. See Vidyasagar |18| equations 7.4.20-21. 

□ 


Proposition 13. The SFB linearizing transformation JT,a,j3 of a SISO de¬ 
terministic dynamical system Qd = if^d^U^xo) G Xn{n,l) into a control¬ 
lable linear system exists if and only if there is a solution A: M"" —> R to the 
set of partial differential equations: 


{d\, adj g) = 0 for 0 < i < n — 2 

{dX,a<i^~^g) 0 . 
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( 21 ) 

( 22 ) 






The linearizing transformation T{x) is given hy: 


■i 

II 

for 

1 < i < n 

(23) 

-CfX 

=- J 1 

C,C}-^X 


1 

C,C}-^X 

(24) 


Proof. See Vidyasagar |]^] equations 7.4.23-33 and Nijmeijer and van der 
Schaft 1^1 Corollary 6.16. □ 


Finally, the geometrical conditions for the existence of the linearizing 
transformation are studied. 


Theorem 14. A deterministic SFB linearizing transformation of Qd = 
{f,g, U,xo) G X£)(n, 1) into a controllable linear system exists if and only if 
the distribution := span {adj 5,0 <i<n — l} is nonsingular on U and 
the distribution A„_i := span | ad j 5,0 <z<n — 2} is involutive on U. 

Proof. See Nijmeijer and van der Schaft |^] Corollary 6.17, Vidyasagar 
Theorem 7.4.16, Isidori Q Theorem 4.2.3 . □ 


2.3 SCT Linearization 


Theorem 15. There is a SCT linearizing transformation Tt of a determin¬ 
istic MIMO system Qd = if, 9 ,U,xo) G X£)(n,m) into a controllable linear 
system if and only if there exists a reordering of the vectorfields 51 ... 5 m 
and an m-tuple of integers < ^2 < ... Km with ^ called the 

controllability indexes such that the following conditions are satisfied for 
all X G U: 

(i) dim ^span |(ad'j. gi{x),Q < i < m, 0 < 5 < k* — 1 )|^ = n 

(25) 

(a) [adj 5 j, adj gj] = 0 for 0 < A: + Z < k* + Kj — 1 , 0 < < m. 

(26) 


Proof. See Nijmeijer and van der Schaft |1^| Theorem 5.3 and Corollary 
5.6. □ 


The following corollary can be verified for SISO systems: 

Corollary 16. For a SISO system with m = 1 the condition (ii) of^ can 
be simplified as follows: 

[ 5 , 045 ] =0, / = 1,3,5, ...,2n-l,Vx G 17. (27) 
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Proof. See Nijmeijer and van der Schaft |13| Corollary 5.6 and the text which 
follows. □ 


3 Transformations of Ito Dynamical Systems 

The transformation rules of Ito systems are motivated by the Ito differential 


rule (see eg. Wong and Hajek |2C] Section 3.3), which defines the influence 
of nonlinear coordinate transformations on Ito stochastic processes. 

The Ito differential rule applies to the situation where a scalar valued 
stochastic process xt defined by a stochastic differential equation dxt = 
f{xt)dt + a{xt)dwt with / : M —> R and ct: R —!■ R smooth real functions 
and wt a Brownian motion, is transformed by a a diffeomorphic coordinate 
transformation T: R —> R. Then the stochastic process zt := T{xt) exists 
and is an Ito process. Further, the process zt is the solution of the stochastic 
differential equation 


dT dT I d^T 

dzt = dt + 


(28) 


All details together with a proof are available for example in Karatzas and 
Shreve |§]. 

The Ito rule can be also derived for the multidimensional case: for the m- 
th component of an n-dimensional stochastic process the Ito rule can be 
expressed as follows: 


2 = 1 


. . I 1 ...in.. 


d^Trr 


2 dxiXj 
*j=i ^=1 


= CfTm dt + dwj + P^Tm dt. (29) 

i=i 

For the most common case with scalar noise k = 1 the equation can be 
further simplified to: 

d'^T 

™ ’■ (30) 


{] z — 


1 "" 

CfTm dt + C„Tm dw + - 


2 ^ 

i,]=i •> 


The operator P^Tm is sometimes written using matrix notation as: 


PaTm = -trace ( a 


d‘^T„ 


2 ^ dx"^ 

Generally, P^ vanishes for linear T or zero a. 


(31) 
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3.1 The Correcting Term 

In this section we introduce an extremely useful equivalence between Ito 
and Stratonovich systems, which allows to use some Stratonovich lineariza¬ 
tion techniques for ltd problems. The motivation is following: let 0/ = 

(/, g, a, U,xo) G X/(n, m, k) be an Ito system. We are looking for a Stratonovich 
system Qs = such that the trajectories of 0/ and 0^ are 

identical. The aim is to find equations relating the quantities /, g, and a 
with /, g, and a. 

Definition 17. Let 0i/ = (/, g, a, U, xq) G X/(n, m, k) be an n-dimensional 
Ito dynamical system with /c-dimensional Brownian motion w. The vector- 
field corro-(x) whose r-th coordinate is equal to 


(corro-(x))r 


1 n ^ Ci 

i V -TV "T (J(Jrj 

Z=lj=l 


for 1 < r < n 


(32) 


is called the correcting term. Note that the derivative is always eval¬ 
uated in the corresponding coordinate system. Further, let us define the 
correcting mapping Corr^ : X/(re, m, fc) ^ X 5 (n, m, fc) by 


CoTv^if, g, cj, [/, xo) := (/ -f corr^(x), c?, cj, U, xq). 


(33) 


The general treatment of the subject can be found for example in Wong 
and Hajek p. 160 or in Sage and Melsa |17|. The following theorem 
describes the behavior of the correcting term under the coordinate transfor¬ 
mation. 


Theorem 18. Let 0/ = {f,g,a,U,xo) G X/(n, 1,1) be a one-dimensional 
ltd dynamical system. Let T be a diffeomorphism defined on U and the sym¬ 
bols Trj, and Trp denote a ltd coordinate transformation and a Stratonovich 
coordinate transformation induced by the same diffeomorphism T and a = 
T*cr. Then the following diagram commutes: 

027 

Corr-fi (34) 

02S 


rf 

01 


ni 


Corrc 


Qis 
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In other words: 


(35) 

(36) 


Trf = Corr^ ^ o o Corio- and 
Trp = Corro-”^ oT^ o Corr^ . 

The notation Corr“^ is used to denote the inverse mapping 
Coj:r^~^{f,g,a,U,xo) := (/- corr^(x), g, cr, C/, xq). (37) 


Proof. The correcting term corro-(x): M”" ^ M"" is equal to 


corro-(x) 


1 da 

2 dx^ 


(38) 


The transformations specified in the diagram (|^ ) will be evaluated in the 
following order: 


017 - > 027 

(b) 

015 - - > 025 > 037 


(39) 


We want to prove the equivalence of 027 and 037. The symbols (a) and 
(c) denote Ito coordinate transformations; the symbols (b) (d) stand for the 
correcting mapping and its inverse. Note that the systems 027, 025, and 
037 are defined in the z-coordinate systems. Further, let 


dT 

T ( ^ dT Ida 
K := T* corr^(x) = 

dx 2 dx 

zo ■■= T(xo). 


(40) 

(41) 

(42) 


(a) 

(b) 

(c) 

(d) 


027 = 
015 = 
025 = 
037 = 


fdT ^ 1 dT 

, cr,U , Zo 

\dx 2 dx^ dx 

r i da 

f - g, cr,U , xo 

2 dx 


dT 


f _ ^ da 

dx \ 2 dx^) ' dx^ ' dx 


dT dT 

g , -T-a ,U,zo 


dT ^ Ida dT 

’ ir9 ,(^^U,zo 

dx 2 dz dx 


(43) 

(44) 

(45) 

(46) 
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All the terms in (a) are equivalent to the respective terms in (d) ex¬ 
cept for the drift terms containing functions of a. Therefore, we continue 
comparing these terms only. For (a): 



92 r 

dx'^ 


For (d): 


(47) 


„ Ida^ Id (dT \ dT Idx d (dT \ dT 


K + 


1 


92 T dT da 
dx‘^ ^ ^ dx dx 


la^T 2 1 2 


Thus L = R and © 2 / = © 3 /. 


(48) 


□ 


Theorem is valid also for combined transformations: 

Corollary 19. Let ©/ = {f,g,a,U,xo) G X/(n, 1,1) , T, T^, and have 
the same meaning as in Theorem^^. Then the following diagram commutes 
for arbitrary regular feedback iFa,g: 


© 1 / -^ ©27 - ^ © 4 / 


CorrcT 


Corr^ 


Corr^^ 


(49) 


7^ Tvv 3 

©IS -*• ©2S -^ ©45 


Proof. We want to prove equivalence of © 4 / and Corr^^^^© 45 . 

The control and dispersion vectorfields of ©47 and ©45 are identical and 
they are not influenced by the correcting mapping. 


Using the notation of Theorem |1^ we can express the drift term of ©47 
as T*/ + L + ga. The drift term of CorryJ ^^©45 is T^f -|- A -|- ga. The effect 
of feedback is purely additive and both the systems are equal. □ 

At first glance the correcting term is rather surprising. How can the sec¬ 


ond derivative of T in (^) be compensated by the correcting term, which 
does not contain the T at all? The answer is quite simple: the second deriva¬ 
tive is hidden in the correcting term implicitly because the correcting term 
depends on the coordinate system in which the system ©17 is dehned. The 
derivative contained in the correcting term is always taken in the appro¬ 
priate coordinate system. To emphasize the dependence of the correcting 
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term on the coordinate system, we will never omit the independent variable 
(eg. X or z) from the symbol corrCT(x). 

Proposition]^ is valid for general multidimensional systems 0/ = (/, g, a, [/, xq) € 
\i{n,m,k) ; the proof is purely mechanical and is not presented here. 

Let us now turn our attention to several special cases of the correcting 
mapping. 

Corollary 20. Let Qj = (f,g,a,U,XQ) G Xj-(n, m, 1) be an n-dimensional 
stochastic dynamical system with an one-dimensional Brownian motion w. 

The r-th coordinate of the correcting term (corro-(x))r is equal to 


(corrCT(x))r. 



i=l 


dUr 1 

dxi 2 


Proof. Substitute k 


1 into (^). 


for 1 < r < n. 


(50) 

□ 


Corollary 21. For systems with one-dimensional noise (k = 1) define the 
matrix valued ltd term PuT for T: M"" —> M"' with components Ti, 1 < i < n, 
as a column n x 1 matrix P^T := [P„Ti,PuT 2 ,..., 

Then the relations ( p0[ ) can he expressed as 

P„T = T* (corro-(x)) - corr 5 -(z). (51) 

Proof. The proof is almost identical to that of the multidimensional variant 
of Corollary (|l 8 |).The symbols can be identified as follows: 

Li = {PaT)i for I <i <n (52) 

Ki = (T* (corr^(x)))j (53) 

Ri = (T* (corr^(x))) • - (corr£^( 2 :)) ■. (54) 

□ 


Corollary 22. Assume that the conditions of Proposition M hold. The 
relation ( 0 ) can he written as: 

]-C^Ca^Ti = PcrTi - Ccorr„{x)Ti for I < i < n. 


Proof. The formula can be expressed as: 


-CaC^Ti = 


2 \ ^dxi ^ I 2^ 

\j=i ^ ) 1=1 


d 


dTi 


0 2-^ ^ dxk * 2-^flnr-, 


— dxj 


k,j=l 


(7j (7 k 




j , dTi daj 


+ 


dxkdxj dxj dxk 


(55) 


CTfc = PaT - £corr,(x)7i- (56) 
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□ 


3.2 Composition of Coordinate Transformations of Ito Sys¬ 
tems 

The set of all deterministic coordinate transformations Tt together with 
composition T^s '■= forms a group. Obviously, this fact is a straight¬ 

forward result of the behavior of the contravariant transformation and there¬ 
fore an analogous statement must hold for Stratonovich systems. Surpris¬ 
ingly, this is valid also for Ito systems as will be shown here. This has an 
important consequence; we may always find the inverse transformation to a 
given coordinate transformation of Ito systems. We will prove the following 
assertion: 

Theorem 23. Let T^, Tg G X/(n, 1,1) be coordinate transformations of 
one-dimensional ltd systems induced by diffeomorphisms R and S. Then 

TioTi = T5V (57) 

Proof. We will transform the system in two different ways and show that 
the results are equal. 

1. In the first method the system A = (0,0,a,U,xo) which corresponds 
to a differential equation dx = a{x) dw will be transformed twice; 

(a) first, by y = R{x) to y coordinates 

(b) and then the result B = {g, 0, b, U, xq) which corresponds to dy = 
g{y) dt b{y) dw by z = S{y) to z coordinates. 

2. The other method transforms the system A only once by z = T{x) = 
S{R{x)) = {SoR){x). 

Without loss of generality, the equation dx = a{x) dw is assumed to have no 
drift term because the drift term transforms in the contravariant fashion. 
The derivatives will be denoted by =; T', =; R\ =: S' and 

similarly for T", R" and S". Note that the prime is always used to denote 
derivatives by the argument of the function. 

The transformation by R gives: 

dx = a dw 

dy = R'a dw -b dt. 
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(58) 

(59) 



Thus the coefficients of the second SDE are defined by 
b{y) := {R'a)oR-^{y) 

9{y) ■= {]^a^R!')oR-^{y). 

The second transformation (by S) gives 
dz = g + S'bdw = 

= + ^{R'fa'^S''^ dt + S'R!adw = 

= (S'R" + (R'fS") dt + S'R'adw = 

The last equality follows from the fact that T" = S'^R')"^ + S'R”. 
One can verify the multidimensional case in the same spirit. 


( 60 ) 

(61) 

(62) 

(63) 

(64) 

(65) 

□ 


3.3 Invariants 

In the deterministic case, some useful propositions about the invariant prop¬ 
erties for example the Leibniz rule (^) and the relation ([l^ were employed. 

Unfortunately, we have not found any analogy for the Ito systems yet. 
To point out the main complications, we will analyze the Ito equivalent of 
the Leibniz rule (0), which is essential for reducing the order of partial 
differential equations in the deterministic exact linearization. 

If the Lie derivative Cg is interpreted as a general first order operator 


4 = E 



( 66 ) 


then the commutator of two such first order operators CfCg — CgCf is also 
a first order operator (see (|l^)). 

Similarly, define the general second order operator as 


0{9,G) 


dxidxj 

i=l * i,j=0 * ^ 


(67) 
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where gi, Gij: M” —> M for 1 < i,j < n. We can compute the commutator 

C(/, F, g, G) := 0(/, F)0{g, G) - 0{g, G)0{f, F) (68) 

of such second order operators. If this commutator was also a second order 
operator (i.e. there were (p and <I> such that the operator G{f,F,g,G) = 
then we would be able to simplify any PDEs of stochastic trans¬ 
formations. (See Proposition 

Because the operator O is linear, i.e. 0{f,F) = 0(/, 0) -|- 0(0, E), we 
can split the computation into four independent, reusable parts: 


Gif, F, g, G) = 0(/, F)Oig, G) - Gig, G)0(/, F) = 

= + 0(0, F, 0, G) + 0(/, 0,0, G) + G(0, F, g, 0). (69) 

The first term is already a first order operator. Only the second and the 
third terms need to be computed because the fourth term can be obtained 
from the third one by formal substitution. For the third term: 


2=1 


o(f,o)0(o,a) = 


92 


dxkdxi 




2,fc,/ = l 


^3 


dxi dxkdxi 


dxidxkdxi J 


Further, 


0(0,G)0(/,0) = Y^Gki 

k,l=l 




dxkdx 


Ef‘ 


A 

dxi 


E o-' 




2=1 

d^h d ^ dfi 52 


+ 2 


+ fi 


33 


dxkdxi dxi dxi dxidxk dxidxkdxi 


(70) 


(71) 


The intermediate results for the third and the fourth terms can be com¬ 
bined into 


0(/, 0)0(0, G) + 0(0, F)Oig, 0) - Gig, 0)0(0, F) - 0(0, G)0(/, 0) = 
dGki dFki . dgi dfi\ d‘^ 


E 

2,Ai,^ = 1 


h 


dxi 


dFki dgi dfi 

9i T ki n 

OXj OXj OXj 


dxkdxi 



dxkdxi 


dxkdxij dxi' 


(72) 
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All of them are first and second order operators. Now let’s evaluate the 
second term 


rP- 

0{Q,F)0{Q,G)=Y,F^3 


dxidxj 




2,J',/C,Z = 1 


^ q2 

^ dxkdxi 

dGki 


dxidxj dxkdxi 


dxj dxidxkdxi 


+FijGki 




dxidxjdxkdxi 


( 73 ) 


0(0, F)O{0, G)-O(0,0)0(0, F)= 


ij,k,l=l 


^ d^Gki ^ d^Fki \ 

dxidxj dxidxj ) dx^dxi 


+ 






dxidxkdxi 


(74) 


Unfortunately the last term 




^3 


dxj J dxidxkdxi 


(75) 


is of third order and, in general, it does not vanish. Thus we have shown 
that the commutator of two general second order operators is of third order. 
Consequently, the Leibniz rule simplifications used in the deterministic case 
cannot be applied to the general stochastic linearization problem. 


4 Stochastic Case 

Since there are two definitions of coordinate transformation of stochastic 
differential equations (Ito , Stratonovich) and three definitions of linearity 
{g,a,ga), we face at least six stochastic problems per a deterministic one. 
In this section we will discuss all of them giving at least partial solutions to 
the feedback linearization problem. We consider mainly the SISO problem 
except for cases where the MIMO extension is trivial. 

4.1 Stratonovich ^f-linearization 

We show that the method known for deterministic systems can be applied 
without modifications. 
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Figure 3: Asymmetry of SFB ^cr-Linearization 

Proposition 24. The Stratonovich dynamical system Qs = if, 9, o") U, xg) G 
Xs(n, m, k) is g-linearizable if and only if the deterministie system 0d = 
{f,g,U,xo) G X£)(n,m) is linearizable. These two linearizing transforma¬ 
tions are equal. This holds both for SISO and MIMO systems and both for 
SFB and SCT linearization. 

Proof. The comparison of transformation laws for deterministic and Stratonovich 
systems shows that the coefficients of / and g transform in the same way. 

The controllability conditions and the definition of linearity are also identical 
(compare Definition ^ with Definition . Identical problems have identical 
solutions. □ 

4.2 Stratonovich ^rcx-linearization 

The Stratonovich problems are not complicated by the second order Ito term. 

The transformation laws for Stratonovich systems are the same as the de¬ 
terministic transformation laws, therefore many results of the deterministic 
linearization theory can be used. 

For example, the stochastic SCT 5 fcr-linearization of a Stratonovich sys¬ 
tem is equivalent to the linearization of a deterministic, non-square MIMO 
system with two inputs and a single output. The SFB problem, which is 
studied in this section, is not as simple as the SCT one because the feedback 
influences only the control input u (Figure ^). The “dispersion input” is not 
a part of the feedback. Consequently, in order to solve the Stratonovich SFB 
5 cr-linearization, we have to deal with a combined deterministic SFB-SCT 
problem. 


4.2.1 Canonical Form 

Recall that we require gf-controllability of the resulting system 
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Since this is a Stratonovich problem, the transformed vectorfields / and g 
do not depend on the dispersion vectorfield a. Therefore, the control part 
and the dispersion part can be studied independently. 

Any g-linear system can be transformed into integrator chain by a combi¬ 
nation of a linear coordinate transformation and linear feedback. Therefore, 
if we set a = 0 , the canonical form is the integrator chain. 

In general, the dispersion vectorfield a is assumed to be arbitrary con¬ 
stant vectorfield d'{x)i = Si, 1 < i < n (See Definition ^ and this form is 
preserved by arbitrary linear transformations. Therefore the canonical form 
can be written as: 


fi{x) ■- 

— Xi-f-l 

for 

1 

< i 

< n — 1 

(76) 

fn{x) ■- 

= 0 





(77) 

9 i{x) ■- 

= 0 

for 

1 

< i 

< n — 1 

(78) 

9n{x) ■- 

= 1 





(79) 

d-i(x) ■- 

= Si 

for 

1 

< i 

< n. 

(80) 


We can compare this canonical form with the equations which dehne the 
transformed system 0. 


Proposition 25. There is a SFB ga-linearizing transformation of the SISO 
Stratonovich system Qg = {f,g,(^,U,xo) G X 5 (n, 1,1) into a g-controllable 
linear system if and only if there is a solution A: M"’ —> R o/ the set of partial 
differential equations: 


{dX, adj 9 ) = 0 

for 0 < i < n — 2 

(81) 

{dX, ad”"^ 9 } 7^0 


(82) 

{dX, ad} a) = s'+i 

for 0 < i < n — 1 

(83) 


sueh that s' G R are eonstants on U for 1 < i < n. Then the linearizing 
transformation is given by: 


Ti = 

-CIX 

CgCy^x 


for 1 < i < n 

( 3 = - 

CgCy^x 


(84) 

(85) 


Proof. Assume that Qs is transformed by JT,a,g into 0 := 


(/,5,d,r(c/),r(xo)) 


where the i-th components of f,g, and a can be expressed as: fi = T/Ti, 
cji = CgTi, di = CaTi- Moreover, the feedback is defined by u = a + jdv. The 
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equations of 0 can be compared to the equation of the canonical form (^)- 


CfTi = Tj+i for 1 < i < n — 1 (86) 

CgTi = 0 for l<i<n — 1 (87) 

CgTn = 1//3 ^ 0 (88) 

CfTn = -a//3. (89) 


The relations (0), (0), (|^), and (0) are equivalent to relations (^)- 
( p^ from Proposition m- The relation (|8^ can be verified in a similar 
way: 


CaTi = Si for 1 < i < n 


thus by (p6| ) 


CfjCfTi = Sj+i for 1 < i < n — 1 


and by 


dD 


CaCfTi = CfCaTi - for 1 < i < n - 1 

since the Lie derivative of a constant is zero: 


(90) 


(91) 


(92) 


CfCcrTi = CfSi = 0 (93) 

Si+i := CaCfTi = for 1 < i < n - 1. (94) 

The equations (|^) are obtained by successive application of this relation. 
The symbols Si are equal to except for the signs. □ 


4.2.2 Conditions for the Control Part 

The necessary conditions for linearizability of the control part of 0 (i.e. the 
system {f,g,0,U,xo)) can be expressed in geometrical form. We intention¬ 
ally omit the dispersion part, using the fact that the resulting system must 
be linear when the noise is zero. 

Further, the class of all solutions of this subproblem will be called C. 
This class can be studied to find if some member of C linearizes the disper¬ 
sion part of the system. 

We can find a geometrical criterion similar to the conditions of Proposi¬ 
tion O. In this case these conditions are necessary but not sufficient since 
also (|83D must be satisfied. 


25 





Proposition 26. SFB ga-linearizing transformation of the Stratonovieh 
system ©5 into a ga-eontrollable system linear system exists only if the dis¬ 
tribution {adj 5 , 0 < i < n — 2 } is involutive and the distribution {adj < 7 , 0 <z<n — l} 
is n-dimensional. 


Proof. This theorem is equivalent to Proposition |T^ which is a direct conse¬ 
quence of Proposition 12 which corresponds to Proposition ^ □ 


4.2.3 Condition for the Dispersion Part 

The conditions of Proposition ^ can be written in matrix form. We are 
looking for Ti = A: M"’ ^ M such that 

adjg 
adjg 

_ ad""^ g 

The vectors adj g... ad^”^ g are written in coordinates as 1 x n rows. 
The first matrix is (n — 1) x n. Moreover it is required that 

{dX,ad'J~^g) (96) 

is nonzero. 

We will use the algorithm for SFB deterministic linearization (see Sec¬ 
tion to find such a transformation A.Then we will verify if the conditions 
for linearity of the dispersion part of the system (^) are also valid. There 
are n additional linearity conditions (sj are constants): 


ax 

dxi 

ax 

ax2 


ax 

axn 


= [ 0 ]. 


(95) 


adj a 


r ax -1 
axi 


•Sl 

adj a 


ax 

ax2 

= 

•52 

adj“^ a 


1 

1_ 


_ _ 


(97) 


In the deterministic case we were satisfied with arbitrary solution A to 
^ 1 ) 


the equations (95), and (|^) . In this stochastic case we must find the class 
of all solutions and then check if this class contains the solution for the a 
part (97). Details depend on the methods used for solving the set of PDFs. 
This result is summarized in the following algorithm; 


• Find Afc := adj s' for 0 < i < A: — 1. 
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• Verify that dim(A„) is n. 

• Verify that A„_i is involutive (see Nijmeijer and van der Schaft [^] 
Remark following Definition 2.39) otherwise no linearizing transfor¬ 
mation exists. 


Find all A satisfying (| 100 |) by solving PDFs (|100D ; denote C the set of 
all such functions. 


• Verify that there is a Ai G C such that the conditions (|^ are satisfied, 
otherwise no linearizing transformation exists. 


Compute T,a,f3 from (p 6 |) - (p9[). 


Now, we can illustrate one possible practical approach which worked for 
several simple problems solved by us (see the example in Section ^). 

First we can compute the kernel of the matrix Mg to find the form uj = 
[ijJi,ui 2 , ■ ■ ■ which satishes MgOJ = 0, i.e. oj is perpendicular to Mg. In 

modern computer algebra systems there is a single command for this. 

Proposition ^ assumes that n vectorhelds A„ := {adj 5 , 0<i<n — 1} 
form an n dimensional space. The vectorhelds A„_i := {adj 5 , 0 < f < n — 2} 
are chosen from them and consequently must form an (n — l)-dimensional 
space. Thus their kernel d\ is exactly one dimensional and arbitrary uj' = 
c{x)u}{x) also belongs to the kernel (c(x) is a scalar). 

But not every oj' that is perpendicular to Mg is a solution to the original 
linearization problem. The function u)' must be an exact one-form i.e. there 
must be a scalar function A such that d\ = c{x)uj{x). The Frobenius theorem 
guarantees that if A„_i is involutive, then there is always c(x) G M such 
that c{x)uj{x) is the exact one-form. 

A necessary condition for a one-form lo = ^o be exact is 


doJi dujj 
dxj dxi 


for 1 < < n. 


(98) 


Hence for every 1 < i,j <n 


{c{x)uJi) = {c{x)ujj ), 


dxj 


(99) 


thus for every 1 < i, j < n 


dc{x) 


dc{x) 

dxj 


LOi + C(x) 




( 100 ) 
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The later condition is a set of linear PDEs, with unknown c{x), which are 
guaranteed to have a solution by the involutivness of A„_i (the Probenius 
theorem). 

In our computations the equation o was in a simple form which 
allowed to determine all the solutions easily. More complicated cases will 
require more sophisticated analysis. 

4.3 Ito g'CT-linearization 

In the previous subsection we tried to find gicr-linearizations for Stratonovich 
dynamical systems. Once this is done, the correcting mapping can be used 
to construct Ito gicj-linearizing transformation. This method works for both 
the SFB and the SCT case. 

Given an Ito system 0/, the corresponding Stratonovich system ©5 can 
be obtained using the correcting mapping ©5 = Corro-(©/). Afterward, 
the Stratonovich gor-linearization algorithm can be applied giving a linear 
system © 25 - Due to linearity of the drift vectorfield a of © 2 s, the correcting 
term coTTa{z) of the backward transformation Corr^-”^ vanishes. 

Theorem 27. The SFB ga-linearizing transformation Ji of the ltd dynam¬ 
ical system ©7 = {f,g,a,U,xo) G X/(n,m,/c) , /(xq) = 0, corro-(xo) = 0, 
into a ga-controllable linear system exists if and only if there is a SFB ga- 
linearizing transformation Js of the Stratonovich dynamical system 


©5 = (^f,g,cr,U,Xo^ = Corr<^(© 7 ) ( 101 ) 

/= / + corr^(x). ( 102 ) 

Moreover Ji = Js o Corro-. 

Proof (sufficiency). We use the properties of the correcting term (Subsection 
0 )- Assume that there is a mapping Js which transform ©5 into a linear 
^-controllable system {Ax,B,S,U,Q) . By (^5|) 

Ji = CorrG^ 0 J 5 o Corr^ . (103) 

The backward correcting transformation CorrG^ is identity because the cor¬ 
recting term of a linear mapping corr 5 -(x) is zero. Thus Ji = Js o Corro- 
and Ji{&i) equals (Ax, B, S,U,0), which is linear and g(-controllable by 
assumption. □ 
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Proof (necessity). Assume that there is the Ito transformation ffi which 
linearizes 0/ and by (|10lD 0/ = Corro- ^( 05 ). Construct Stratonovich 
linearization by Js = JioCoxia~^■ Hence Ji linearizes Corro-“^( 05 ) and Js 
linearizes 0 ^ into the same linear and controllable system as Jj. □ 

4.4 Ito ^'-linearization 

The Ito gf-linearization problem is probably the most complicated variant of 
exact linearization studied in this paper. The dispersion vectorfield of an Ito 
dynamical system transformed by a coordinate transformation Tp consists of 
two terms: the transformed vectorfield T*/ and the Ito term Pq-. We require 
that the sum of these terms is linear, thus the nonlinearity of the drift T*/ 
must compensate for the Ito term. Since the Ito term behaves to T as a 
second order differential operator, this problems generates a set of second 
order partial differential equations. One can attempt to use simplifications 
as in the deterministic linearization, namely, the recursive Leibniz rule ([l^. 
Unfortunately, this approach does not work for the stochastic case. In gen¬ 
eral, the Ito equations cannot be easily simplified by commutators because 
the commutator of second order operators is not a second order operator 
but a third order operator (see Subsection |3.3|) . 

Nevertheless, there are special cases for which simpler conditions can 
be found. The most important special case (commuting g and cr) will be 
studied here. 

4.4.1 Canonical Form —n unknowns 

The canonical form for the ^r-linearization is the integrator chain with a 
nonlinear drift 


fi{x) = Xi+l 

for 1 < i < n — 1 

(104) 

fn{x) = 0 


(105) 

m{x) = 0 

for 1 < i < n — 1 

(106) 

gn{x) = 1 . 


(107) 


Assume that there is a ( 7 -linear system 0/ = {Ax, B,a{x),U,xo). Then 
the drift part of 0 / can be transformed by a linear transformation into the 
integrator chain. This is because the Ito term of a linear transformation 
vanishes. 

The equations which define T can be obtained by comparing this canon¬ 
ical form with the equations of 0 . 
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Proposition 28. Let 0/ = (/, g, a, U, xq) £ X/(n, 1,1) be an ltd dynamical 
system with /(xq) = 0 such that corro-(xo) = 0. There is a SFB g-linearizing 
transformation JT,a,i3 of the system 0/ into a g-controllable linear system 
if and only if there is a solution Ti : M” —> M, 1 < z < n, to the set of partial 
differential equations defined on U: 


Tj+i = CfTi + PfjTi for I <i <n-l (108) 

CgTi = 0 for 1 < i < n — 1 (109) 

CgTn / 0 . ( 110 ) 


The symbol denotes the ltd operator (see Definition^). The feedback can 
be constructed as: 


a = — 


{k^fTn + PaTn) 


C T 


/9 = 


T T 


( 111 ) 


■jii -f- gi 


Proof. The i-th components of f,g, and a are: jt 
di = CfjTi. The partial differential equations 
obtained by comparison of (^)-(^) with the equations (104)-(107). 




(|10g|) and (|110D are 
□ 


4.4.2 PDEs of single unknown 

0ne can attempt to reduce the equations (108), (109), and ( |110D , to a set 
of equations of a single unknown, similarly to the results of Proposition 


Corollary 29. Define the general second order operator 0{f,F) as in ( |t 
The exponential notation for O will be defined recursively: 0^{f,F)T := T 
and 0^+^{f,F)T := 0{f,F)0\f,F)T for I > 0. Next, define 


F - — 


( 112 ) 


Then the set of partial differential equations (|108|) - (|110D ofn unknowns has 
a solution if and only if there is a solution A: M”" —> M defined on U to the 
set of PDEs of single unknown: 


0{g,0)0\f,F)X = 0 for 0<i<n-2 (113) 

O{g,0)O^-\f,F)X^0. (114) 

The original solution and the feedback can be found as 

Ti = 0^-\f,F)X for l<i<n 

0{f,Fr-^X 1 

O(5,0)-iA ^ O(5,0)-iA' 
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(115) 

(116) 














(117) 

(118) 


Proof. Since Ti = A and by definition of (^) and O (^: 

0{f, F)Ti = CfTi + P„Ti for 1 <i <n 
0{g,0)Ti = CgTi for 1 < i < n; 

then Ti+i = 0{f,F)T, by (|T^) and T, = 0^-\f,F)n = 0^-^{f,F)X. 

Similarly the equation ( |111|) implies Proposition |116| . □ 

Note, that for the deterministic case a = 0, the operators 0{g, 0) and 0{g, 0)0® (/, F) 
degenerate to Cg and to FgD^ respectively, thus the result is the same as 
that of Proposition 

In general, the equations of the system are of an order up to 2n and 
cannot be reduced to a lower order. The commutator of two second order 
operators is of third order as was pointed out by (^5|) . In particular, for i = 1 
we have to evaluate C{g, 0, /, F) =: 0{a, A), which is of second order due to 
the fact that G = 0. But starting from i = 2 the commutator C{f, F, a, A) = 

C{f, F, C{g, 0, /, F)) is of third order. 

4.4.3 Correcting Term 

The same problem can be reformulated using conversion to the Stratonovich 
formalism. One can compute the drift vectorfield / of the equivalent Stratonovich 
system by applying the correcting term f ■= f + corro-(x). Then the 
Stratonovich system (^f, g, a, U, may be transformed, by a suitable trans¬ 
formation, to such a form (/, g, d, T(t/), 0) that after applying the backward 
correcting term — coviaiz) the resulting Ito system will be linear. 

Compare this formulation with the ^c-linearization where the backward 
correcting term vanished due to linearity of a. This does not happen with 
the ^(-linearization, and the backward correcting term is a part of the equa¬ 
tions. 

In general, it may be difficult to solve these equations. Nevertheless 
there are special cases when the solution can be obtained. See for example 
Section ||. Another important case (commuting g and a) is studied below. 

Corollary 30. The equations of Proposition Wh are equivalent to n partial 
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differential equations: 


Tj+i = CjTi- coTTs{z) = (119) 

C^Ti + for 1 <i <n-l (120) 

CgTi = 0 for l<i<n-l (121) 

CgTr, / 0. (122) 


Then the feedback can be constructed as 


C -T ^ -C C T 

= - - -;r7^;- 

C 1^ 

n 


/3 = 


r T 


(123) 


Where / := / + corro-(x). 


Proof. The equations 
and 


can be obtained from ( |108D by applying (|^ 


Ti+i — CfTi + P„Ti — + Po-Tj — 

— Ccorrg{x)'^i + Pa^i = P fPi + -^PaPaPi- ( 124 ) 

The other equations are adopted from Proposition The set of PDEs 
of n unknowns can be transformed into a set of PDEs of a single un¬ 
known A = Ti, but the order of the equation will be 2n — 1. □ 


Remark 31. Observe that the set of n second order partial differential equa¬ 
tions (119)- (122) defined in Proposition ^ can be transformed, by intro¬ 
ducing new variables Si = C^Ti, to the following system of 2n — 1 first order 
partial differential equations for 1 < i < n — 1: 


CgTi = 0 (125) 

CaTi - Si = 0 (126) 

Cfi + ^C^Si = Ti+i (127) 

CgT^ + 0. (128) 

Ti and Si : M”" —> M are unknown real valued functions defined on t/. V 
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4.4.4 Systems with Commuting g and a 

There is a special case of Ito dynamical systems for which the solution is 
completely known and can be computed using only first order PDEs. 

Theorem 32. Let 0/ = (/, g, a, U, xq) G X/(n, 1,1) be a SISO ltd dynamical 
system. If the vectorfield a commutes with all vectorfields adLg for 0 < z < 

n — 1, i.e. [adyg,a] = 0 , where / := / + corro-(x), then the ltd system is 
g-linearizable if and only if the distribution 

A„ := span|adEgr,0 < z < n — l| (129) 

is nonsingular on U and the distribution 

A„_i := span |adEg(, 0 < z < n — 2| (130) 

is involutive on U. If these conditions hold, then a solution A: M” —> M to 
the set of partial differential equations exists 


{dX, ad-Lg) = 0 for 0 < z < rz — 2 (131) 

{dX,ady^g) ^0. (132) 

the linearizing transformation is given by: 


Ti = X 


1 


7)+! — d- 

-c^x 

a =- 1 — 

LqL^-^X 

® / 


for 1 < z < n — 1 


P = 


C„CL-^X 


(133) 

(134) 

(135) 


Proof. We will apply the Leibniz rule to the relation of (|119D Corollary ^ 
to expand the term CgTi+i for 1 < z < n — 1 


h^gTi+l — Cg^CjTi + —CfjCaT^ — CgCjTi + —CgCaCaTi — 
C-jC-gTi — + —CaCg{LaTi) — —C[a-,g]{dIaTi) = 


1 


1 


^\Kg]^id- Ca{hlahlgTi Tj^gjTj) (CrjC[a,g]Ti — 


- C 


d- ^L[a,[a,g]]Ti. (136) 


33 






Figure 4: Ito and Stratonovich cr-linearization 


If the vectorfields g and a commute, then the second and third terms vanish. 
If, moreover, [a, [/,(?]] = 0 then 


C-gTi+2 — —C, 


{f,9\ 


CjTi + j - 


(137) 


In general if [a, sdyg] 

CgTk = 


= 0 for 0 < i < n — 1 then 


Thus the equations (|119|) and (|121|) will be equivalent to 


(138) 


{dX, adhg) = 0 for 0<z<n — 2 (139) 

{dX, ady^ g) ^ 0, (140) 


which are of the same form as the equations of Proposition and conse¬ 
quently the conditions from Proposition Id can be used. □ 


4.5 Ito and Stratonovich a-linearization 

The stochastic SFB cr-linearization problem is similar to deterministic SCT 
linearization. The dispersion vectorfield a transforms in the same way as 
deterministic drift vectorfields do. Consequently, no ltd term complicates 
the transformation. Moreover, the Ito and Stratonovich cases are equivalent. 

On the other hand, in the SFB cr-linearization we are free to choose the 
feedback !Fa,i3 that perturbs the drift vectorfield / into f = f + ga. 


Proposition 33. Let Q be a SISO stochastic system 0 = {f,g,a,U,xo). 
There is a SFB a-linearizing transformation JT,a,j3 into a a-controllable 
linear system if and only if there is a smooth feedback function a: M”' ^ M 
such that the deterministic system {f + ga,a,U,xo) has a SCT linearizing 
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transformation Tt- Equivalently, there must be sueh a that the the modified 
odd bracket condition: 

[a,ad^f^g^a] =0 for l = l,...,2n-l (141) 

is satisfied (see (^7|) ). The resulting combined transformation consists of 
composition of the coordinate transformation Tt and the feedback Ta,p where (3 
is arbitrary function of x; for instance (5 = 1. 

Proof. Compare definition of linearity of a deterministic system with defi¬ 
nition of u-linearity. The system is is a-linearizable if and only if the deter¬ 
ministic systems (/ -|- ga, a, U, xq) is SCT linearizable (see Corollary [T^) . 

It is evident that the function (3 (see Figure ^ has no effect on the dis¬ 
persion part and can be chosen arbitrarily. (Probably nonzero for otherwise 
the system will be g(-uncontrollable). □ 

The condition ( |141[ ) can be expressed in terms of derivatives of a using 
bracket relations known from differential geometry. For example, for / = 1: 

[a, [f + ga, cr]] = [cj, [/, a] -h [ga, cr]] = [a, [f, cj]] -h [a, [ga, cr]] = (142) 

= Vr-,g+ [f + g,(T]]+gC„C„a. (143) 

The other conditions for k = 3, 5, 7,... can be expressed in a similar way 
giving the set of n partial differential equations of the order up to 2n for 
example by a computer using symbolic algebra tools. The problem is not 
very interesting from the practical point of view. 

5 Example—Crane 

In this section the methods of stochastic exact linearization are demon¬ 
strated on an example — control of a crane under the influence of random 
disturbances. The description of the plant was adopted from Ackermann 
et al. iJH where the model of a crane linearized by approximative methods 
was studied. Unlike Ackermann, we control the same system using the exact 
model. Moreover the influence of random disturbances is added. 

Consider the crane of Figure which can be used for example for loading 
containers into a ship. The hook must be automatically placed to a given 
position. Feedback control is needed in order to dampen the motion before 
the hook is lowered into the ship. The input signal is the force u that 
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Figure 5: Crane 


accelerates the crab. The crab mass is me, the mass of the load the 
rope length is I, and the gravity acceleration g. 

We assume that the driving motor has no nonlinearities, there is no 
friction or slip, no elasticity of the rope and no damping of the pendulum 
(eg. from air drag). We will define four state variables: the rope angle xi 
(in radian), the angular velocity X 2 = xi, the position of the crab X3, and 

the velocity of the crab X 4 = X3. As shown in ||^], the plant is described by 

two second order differential equations: 

u = {me + mc)x^ + mLl{xi cosxi — x\ sinxi) (144) 

0 = mixz cosxi + m^lxi + meg sinxi. (145) 


Additionally, we assume that the load is under influence of random distur¬ 
bance, which can be modeled as a white noise process. The disturbance 
(wind) is horizontal, has zero mean and can be described by the Ito differ¬ 
ential dw: 


, F cos xi , 

dx 2 = -;— dw, 

mil 


(146) 


where F is a constant having the physical unit of force. 

We used symbolic algebraic system Mathematica to handle the compu¬ 
tations. The complete Mathematica worksheet can be downloaded from the 


web page of the author tittp://www.tenzor.cz/sladecek 


Mathematica was used to solve the equations of the system for unknown 
values X 2 and X 4 (angular and positional acceleration). Values of vector- 
helds f,g, and a were derived as follows: 
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Figure 6 : The State Space Model of Crane 


/ = 
9 = 


sm.xi{g(mL + mc)+ lmLX2Cosxi) 

X2, - 77-^-5 — ^ - ,xa, 0 


n T 


l{mc + rriL — rriL cos^ xi) 


0 ,- 


COS Xl 


l{mc + niL — rriL cos^ xi) 


,0,t6 


Fcosxi 

0 ,-^, 0,0 

rriLl 


n T 


(147) 

(148) 

(149) 


The state space model is shown in Figure |^. We can see that the 
positional state variables X 3 and X 4 are isolated from the angular state 
variables xi and X 2 - Later, we will concentrate on the angular variables 
pretending that the load will be stabilized no matter where the crane is. 
Consequently, we obtain only two-dimensional system on which the exact 
linearization techniques can be demonstrated. 


Next, consider the random disturbances. Because the correcting term corro-(x) 
is zero, there is no difference in using either the ltd or the Stratonovich inte¬ 
gral. In case of more “nonlinear” noise, one of the integrals must be selected. 

If the ltd model is chosen. Theorem ^ must be applied. 

Now we evaluate the conditions of Proposition ^ to check that the 
system is linearizable. In fact, we must only evaluate the non-singularity 
condition because every one-dimensional distribution is involutive, and the 
integrability is satisfied automatically. To this end, we will compute the 
null space (kernel) of the matrix \\f,g\,g\, which is empty and therefore the 
matrix is nonsingular. We conclude, that the deterministic SFB problem is 
solvable. 
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Notice, that the system is already in the integrator chain form and 
hence A = xi satisfies this condition. Therefore, the deterministic system is 
linearizable by feedback only, with no state space transformation at all, i.e. 
z = T{x) = X. 

This choice of the output function A is natural but does not cancel the 
nonlinearity in the dispersion coefficient a. For this purpose, we must use 


the algorithm of Section 4.2.3 to construct another nontrivial coordinate 
transformation T. 

To obtain this transformation, we must find the space of all functions A 
satisfying conditions for feedback linearity (^). Observe that Cg\ must be 
zero hence 


dX dX 

■^91 + -^92 = 0 . 

OXi 0X2 


(150) 


Since 51 = 0 and 52 7 ^ 0 in neighborhood of xq, then = 0 and 
A = ci(xi) is a function of xi only (i.e. without X 2 ). The coordinate 
transformation is T = [A, £jA]^. We want to select such ci(xi) that the 
dispersion vectorfield a := r*iT in the new coordinate system 2 ; = T{x) will 
be constant: 


dci{xi) F cos xi 
dxi mil 


= constant. 


We decided to define the constant as F/{mLl), therefore 

1 


dci 

dxi 


cos Xl 


(151) 


(152) 


and 


Ti — A — Cl 


(xi) = [ 

J COSXl 


■ dxi = 


, / 1 ■ ^ \ ^ ^ -1 
— In cos -Xl — sm -xi + In cos -xi + sm -xi 

'2 2 V 2 2 


(153) 


T 2 = £/A = X 2 secxi. (154) 

Finally, we can compute the feedback from (^5|). In the Mathematica 
worksheet we validate the results by computing 0 = JT,ot,( 3 {f■,9-,(^-,UiXo)-, 
The computation showed that the system ^/, g, a, U, xq^ is in the integrator 
chain form in the 2 ; coordinate chart. 


38 
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I yjic + mi (sm(xi))^ j 

a = tan(a:i) (sec(xi) X 2 — bg {me + mi) + lmiX 2 cos(xi)) . (155) 


6 Conclusion 


6.1 Main Results 

(i) The structure of the stochastic linearization problem is much richer 
than the structure of the deterministic one. Two definitions of coordi¬ 
nate transformation exist and there are differences between the g, a, 
and 5 (T-linearization. 

(ii) In the case of Ito integrals, the coordinate transformation laws are of 
second order (the Ito rule). 

There is a large difference between ^nT-linearization and (^-linearization. 
In the former case the effect of the Ito term can be reduced to the first 
order operator and consequently the problem is solvable by differential 
geometry. On the other hand, in the later case there is no easy method 
to elliminate the Ito term and the a set of second order partial differ¬ 
ential equations must be solved to get the linearizing transformation. 

(iii) We have given (at least partial) solutions to all SISO SFB problems. 
The results are listed in Table |l[ 


Linearization 

Deterministic 

Stratonovich 

Ito 


9 9(y cr 


Theorem 14 

p. 13 









Propos. |24 

p. |23| 

Theorem 

26 

P- 

26 

Propos. 

33 

P- 

34 

Corollary |2f 

^p. ^ 

Theorem 

27 

P- 

28 

Propos. 

33 

P- 

34 


Table 1: Overview of results — SISO SFB case 


(iv) Ito linearization problems can be approached by means of the cor¬ 
recting term. The Ito differential equation can be converted to the 
Stratonovich equation whose behavior under coordinate transforma¬ 
tions is simpler. This method is only partially applicable to the g- 
linearization. 
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(v) An important special case was identified for the ltd ^r-linearization. 
The case is characterized by commuting control vectorfields g and 
dispersion vectorfields a. Solutions can be found using first order 
methods. 

(vi) Computer algebra proved to be a useful tool for solving exact lineariza¬ 
tion problems. 

(vii) Industrial applications of the exact linearization in general are still 
unlikely, mainly due to complexity, sensitivity, and limited robustness 
of the control laws designed by the method. 

6.2 Future Research 

(i) Find a solution to the ltd ^'-linearization problem in general case, 
including geometric criteria, using second order geometry. 

(ii) Analyze the computability issues; implement a universal symbolic al¬ 
gebra toolbox for the problem. 

(hi) Solve the SCT problem. 

(iv) Extend the results to the MIMO systems. 

(v) Extend the results to the input-output problems and linearization of 
autonomous systems. Work out the applications of nonlinear filtering. 

(vi) Perhaps, some of the results can be used as a starting point for ap¬ 
proaching more general class of problems as the problems of distur¬ 
bance decoupling, input invariance of stochastic non-linear systems, 
or problems of reachability and observability. 
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